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Abstract 

In this paper, we review the quantum mechanics of magnetic reso- 
nance imaging (MRI). We traverse its hierarchy of scales from the spin 
and orbital angular momentum of subatomic particles to the ensemble 
magnetization of tissue. And we review a number of modalities used in 
the assessment of acute ischemic stroke and traumatic brain injury. 
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1 Introduction 

MRI plays a key role in the diagnosis and management of acute ischemic stroke 
and traumatic brain injury. A rigorous understanding of MRI begins with the 
elementary particles which give rise to the magnetic properties of tissue. El- 
ementary particles have two independent properties which manifest magnetic 
moments: (i) orbital angular momentum, and (ii) intrinsic spin. These two 
properties are independent only to first order, as it was their very interaction 
via spin-orbit coupling which led to detectable energy differences in the hydro- 
gen spectrum, and consequently to the discovery of intrinsic spin. We discuss 
both properties in this paper in the context of the hydrogen atom. The hydro- 
gen atom is the most frequently targeted nucleus in MRI due largely to its 
biological abundance and high gyromagnetic ratio. The nucleus consists of 
a single proton. Protons are made of quarks, specifically one down and two up 
quarks, each of which are of spin 1/2. The charge and spin of the proton are 
both directly due to its quark composition. The down quark has a charge of -1/3 
while the up quarks each have charge +2/3, which all sums up to a charge of 
-t-1. The nuclear spin is 1/2 by virtue of spin cancellation from the antiparallel 
alignment of two of the three quarks. 




Figure 1: T2 axial fast spin echo image showing acute ischemia in the PC A 
distribution. There is relative hyperintensity in the left occipital region of in- 
farction. Note the prominent hyperintensity of the vitreous and cerebrospinal 
fluid, both consistent with the hyperintensity of water on T2-weighted imaging. 



Content Outline 

The remainder of this paper is organized as follows: Section ^ presents an 
overview of the MRI, Section ([3| reviews the MRI assessment of acute ischemic 
stroke, Section Q describes the hydrogen atom and reviews an analytical solu- 
tion to its Schrodinger equation. Section ([5| reviews intrinsic spin. Section ([6| 
reviews quantum mechanical addition of orbital angular momentum and spin. 
Section ^ reviews Group Theory of SU(2) and SO (3) and their relationship, 
and Section ([s]) summarizes this review. 

The MRI section. Section ([2]), contains a number of subsections which are 
outlined as follows: Subsection ( |2.1[ ) discusses the hierarchy of scales in MRI 
phenomenology. Subsection (|2.2|) reviews the magnetism of electrons and nu- 



cleons. Subsection (2.3) reviews the magnetic potential energy and its role in 



the descriptions of magnetization, torque, force, and larmor precession. Subsec- 



tion (2.4 1 reviews pulse magnetization. Subsection (2.5) describes equilibrium 



magnetization. Subsection (2.6) describes spin-lattice or Ti recovery. Subsec- 



tion (2.7) describes spin-spin or T2 decay. Subsection (2.8) reviews free induc- 



tance decay. Subsection (2.9) reviews the Bloch equations. Subsection (2.10) 



overviews intravenous contrast agents in MRI, Subsection (2.11 ) reviews gradient- 



based band-selection method for MRI space localization, and Subsection (2.12) 
reviews the components of the MRI machine. 

The stroke section, Section ([s]), contains a number of subsections which are 
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outlined as follows: Subsection (3.1 1 reviews diffusion weighted imaging (DWI), 
Subsection (3.2) reviews perfusion weighted imaging (PWI), Subsection (3.3) 
reviews the combination of DWI and PWI and its utility in imaging the penum- 
bra, Subsection (3.4) reviews magnetic resonance spectroscopy, Subsection (3.5) 
reviews blood oxygen level dependent (BOLD) imaging, and Subsection (3.6) 
reviews MRA. 

Appnedix (A) reviews the signals processing involved in MRI, including the 
Fourier transforms, the canonical sampling theorem, and the convolution theo- 
rem. 



2 Magnetic Resonance Imaging 

The MRI is an optical technology whose core equation is the Maxwell-Faraday 
Equation. It exploits ensemble phenomena in which the composition of a sample 
can be probed by sensing its magnetic properties through radio frequency waves. 
Here, we describe the various parts of MRI and their respective roles in the 
whole. 

2.1 A Hierarchy of Scales 

The hierarchy of scales in the modeling of magnetic resonance imaging proceeds 
from the subnuclear scale of quarks and gluons, to the subatomic scale of discrete 
electrons and nucleons, and finally to the bulk sample scale where ensemble 
effects and statistical mechanics are the operative physics. In a general sense, 
the subnuclear scale is coupled to quantum field theory, the subatomic scale is 
coupled to quantum mechanics, and the bulk sample scale is coupled to classical 
electromagnetics. 

Certain phenomena in MRI are only meaningful over specific regimes. For 
instance, a magnetic moment experiences a torque when an external magnetic 
field is applied. However, the concept of torque is a classical concept. It is 
deterministic and continuous, and is meaningful only at the level of ensemble or 
bulk effect. On the microscopic level, such as on the scale of single electrons and 
nucleons, quantum mechanics is the operative physics, and state transitions are 
quantized. This is encoded in the time evolution of quantum states governed 
by an appropriate Hamiltonian, yielding the specific Schrodinger equation that 
applies in that setting. 

Although pedagogically separable, the scales are intrinsically linked physi- 
cally. For instance, the bulk magnetization, M, is the net sum of the atomic 
level magnetic moments, iJby Similarly the pulse magnetization frequency used 
to torque the bulk magnetization vector and change its direction is a radio fre- 
quency wave, Bi, oscillating at the larmor frequency of the individual nucleons. 
Furthermore, Bi takes its effect by acting directly on the individual electrons 
and nucleons, causing them to precess in phase. 

2.2 Magnetism of Electrons and Nucleons 

Electrons, protons, and neutrons are magnets. They derive their magnetism 
from intrinsic spin and orbital angular momentum. In the presence of an ex- 
ternally applied magnetic field, this results in bulk magnetism of sample tissue. 
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The magnetic dipole moment, /x, of a single unpaired electron or nucleon is 
given by, 



11 = 7S, (1) 

where 7 is the gyromagnetic ratio measured in Hz/T (Hertz per Tesla) and 
S is the spin. The gyromagnetic ratio for the electron is given by, 

(2) 



2mp 



where e is the elementary charge, is the mass of the electron, and is the 
electron gyromagnetic factor (g- factor), a dimensionless number whose experi- 
mental agreement with theory is near unprecedented on this scale, and thereby 
lends strong credibility to the theory of quantum electrodynamics. Table ^ 
shows the electron and proton g-factors. The gyromagnetic ratio for nucleons 
is given by, 

where mp is the mass of the proton, and is the g-factor of the nucleon. 
Estimated gyromagnetic ratios of some biologically-relevant nuclei are shown in 
Table Q [13110] ■ 



Table 1: Gyromagnetic ratios of some biologically-relevant nuclei 



Element 


Nuclei 


7„ (10" rad s"' T"') 


7„/27r(MHz T^-^) 


Hydrogen 


^H 


267.513 


42.576 


Deuterium 


D, ^H 


41.065 


6.536 


Helium 


^He 


-203.789 


-32.434 


Lithium 


^Li 


103.962 


16.546 


Carbon 


13q 


67.262 


10.705 


Nitrogen 


14n 


19.331 


3.007 


Nitrogen 


15N 


-27.116 


-4.316 


Oxygen 




-36.264 


-5.772 


Fluorine 


19p 


251.662 


40.053 


Nitrogen 




19.331 


3.007 


Sodium 


2^Na 


70.761 


11.262 


Phosphorus 


31p 


108.291 


17.235 



The Bohr magneton, fiB, and the nuclear magneton, /ijv, are named quan- 
tities related to the electron and proton gyromagnetic ratios respectively, and 
are given by, 

f^B = ^ — , (4) 

and 

The magnetic moment of the electron is much larger than that of the nu- 
cleons, and as is discussed below, this corresponds to smaller precession, pulse, 
and signal frequency for protons than electrons. Specifically, nucleons precess. 
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absorb, and emit cilcictromagnctic signals in the radio frequency range. This 
range is non-ionizing radiation, and makes MRI a safer choice than ionizing 
imaging modaUties such as plane and computed tomography X-rays. 




Figure 2: Tl FLAIR. Sagittal section. Note the hyperintensity of orbital fat, 
and the hypointensity of vitreous and cerebrospinal fluid. Hyperintensity of fat 
and hypointensity of water are radiologic characteristics of Tl weighted imaging. 



2.3 The Magnetic Potential Energy 

The magnetic potential energy is, 

U = -nB. (6) 

The force on a magnetic moment, the torque on a magnetic moment, and 
the larmor precession frequency or resonance frequency, can all be described in 
terms of the magnetic potential energy, U{r,t). 

2.3.1 Magnetization 

The magnetic moment aligns itself in either parallel or anti-parallel orientation 
to the external field Bq. The extent of parallelism is constrained by the eigenval- 
ues of the and operators above, such that from the geometry, the parallel 
and anti-parallel states make angles of cos~^(l/-\/3) « 54.7356103172 degrees 
with the positive and negative z-axes respectively (assuming Bq is oriented ex- 
clusively along the z-axis). 

This process is called magnetization. In the case of a population sample, 
upon application of Bq, the spin orientation distribution goes from a somewhat 
arbitrary array of directions to being composed of only two directions, parallel 
and antiparallcl. The parallel orientation is the lower energy state, and thereby 
is the more populous state for protons in the sample. The population distri- 
bution is temperature dependent, such that the proportion of protons in the 
higher energy state increases with temperature. The relationship is given by, 
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where Nh is the higher energy state, A^; is the lower energy state, E is the 
energy difference between the two states, k is the Boltzmann constant, and T 
is the absolute temperature in Kelvin. The bulk magnetization can be modeled 
as, 

M(t)=^/x^.(0. (8) 
j 

with Bo = zSo, it follows that Af(0) = M^(0) and M^y = 0. In other words, 
the equilibrium magnetization is entirely along the z (longitudinal) direction 
and there is no component in the xy (transverse) direction. 

2.3.2 Torque 

The torque on a magnetic moment in a magnetic field acts to align it with the 
field, is a function of the magnetic potential energy, and is given by, 

r{a) = -^{-,^-B), (9) 

where a is the angle between B and /lx, and the minus sign indicates the 
restoring nature of the torque. 

2.3.3 Force 

The force on a nucleon or electron in a magnetic field is, 

F = -V(-/x-B). (10) 

2.3.4 Precession 

When an external magnetic field Bq is applied to an electron or neutron par- 
ticle, the particle's magnetic moment precesses about the direction of Boat a 
frequency called the larmor frequency given by, 

u; = jBo. (11) 

In the equilibrium magnetization state with Bo = zBq and M^y = 0, there 
is no precession of the bulk magnetization M because /lx = M x Bq = 0. 

Considering only the spin of a spin 1/2 object in a magnetic field, and 
assuming B = e^, we get the following Hamiltonian, 

H = -iJ,-B = -{^^S)-B:=u}-S, (12) 

ZTflx 

where the subscript x denotes either e for electron or p for a neutron, and 
where, 

9xeB^ 

It follows that, 
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H = ujSz. 



(14) 



We see that remarkably, the Hamiltonian is proportional to the z-component 
of spin, therefore the magnetic spin eigenstates are simultaneously energy eigen- 
states. Furthermore the energy eigenvalues are proportional to the magnetic 
spin quantum numbers. Specifically, 

E±^±^huj. (15) 

The above relation led early researchers to the discovery of spin, by its 
manifestation as the Zeeman effect and the splitting of energy levels in the 
hydrogen optical spectrum. 



2.4 Pulse Magnetization 

Given a bulk magnetization M due to a static magnetic field Bq = Bqz, a 
pulse of a weaker magnetic field, Bi, oscillating at the larmor frequency and 
applied perpendicular to Bg, will result in a transition of the individual nucleons 
from dephased to in-phase oscillation. This manifests as the acquisition of a 
transverse phase M.xy and oscillations of M about the z-axis at the larmor 
frequency uj. This "tilting" of the magnetization vector also implies a decrease 
in Mz. During the pulse, the net external magnetic field is the sum of Bq 
and Bi. If the pulse duration was extended indefinitely, the result would be 
alignment of the bulk magnetization with the new direction, Bq + Bi, and there 
would be no resultant oscillation observed. The pulsed nature of Bi is therefore 
critical to the technology of magnetic resonance imaging. 

2.5 Equilibrium Magnetization 

The net magnetization is given by, 

M = (TV, - Nh)tJ,, (16) 

while the ratio of Ni to Nh is given in Equation yn above as, ^ = e~^/'''^ . 
At equilibrium, the solution of the above system oi two equations yields the 
equilibrium magnetization, 

where k = k/2TT is the reduced Boltzmann constant and pp is the proton 
density. 



2.6 Spin-Lattice Relaxation (Tl Relaxation) 

The spin-lattice relaxation, also known as Ti relaxation or longitudinal relax- 
ation, is the process of longitudinal magnetization recovery following a Bi per- 
turbation. Specifically, Ti is the time constant of longitudinal magnetization 
recovery, and is given by, 

Mz{t) = Mo(l - e-*/^i) + M,(0)e-*/", (18) 
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where Mq is equilibrium magnetization and (0) is the longitudinal magne- 
tization instantaneously after pulse excitation. In inversion recovery, the initial 
magnetization is negative, and typically Mz(0) = — Mq. Substituting this value 
into Equation ( 18 ) above yields. 



(19) 



Fluid attenuated inversion recovery (FLAIR) is useful in the diagnosis of 
certain pathologies such as periventricular white-matter lesion in multiple scle- 
rosis. FLAIR can be Tl or T2 weighted. Figure ([2| shows a sagittal section of 
a Tl FLAIR brain image, and Figure ([3| shows a T2 FLAIR brain image. 




Figure 3: T2 FLAIR showing hyperintense signal in the PCA and MCA distri- 
butions, consistent with acute ischemic stroke. Note that the hypointensity of 
cerebrospinal fluid due to fluid-attenuation allows for better contrast with the 
periventricular lesion. 



2.7 Spin-Spin Relaxation (T2 Relaxation) 

The spin-spin relaxation, also known as T2 relaxation or transverse relaxation, is 
the process of transverse magnetization relaxation following a Bi perturbation. 
Specifically, T2 is the time constant of transverse magnetization recovery, and 
is given by, 

M^y{t) = M^y{0)e-'/^^, (20) 
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where Mxy{0) is the transverse magnetization instantaneously after pulse 
excitation. Following the Bi induction of in-phase precession, the proton spins 
again begin to dephase according to T2 time constant. Of note, the relaxation 
of transverse magnetization is faster than the recovery of longitudinal magneti- 
zation. In other terms, T2 <Ti. Figure ([T]) shows a T2 fast spin echo image of 
acute ischemic stroke in the posterior circulation. 



2.8 Free Inductance Decay 

Upon pulse excitation and tilting, the magnetization vector oscillates about the 
z-axis. This oscillating magnetic moment in turn generates an electromotive 
force in accordance with Faraday-Maxwell law of electromagnetic induction, 

OB 

VxE^--, (21) 

where Vx is the curl operator, E(r, t) is the electric field, and B(r, t) is the 
magnetic field. The above equation can be rendered in integral form via Stokes 
theorem to yield, 

E-ds = -^ ff B-dA, (22) 

Ids OiJjT, 

where 9E can represent a wire coil bounding a surface S, ds is an infinitesi- 
mal segment length along the wire, and dA. is vector normal to an infinitesimal 
area of E. Such integral transformation has certain advantages in numerical 
computation. For instance, Green's functions approaches can be invoked and 
can greatly decrease the number of unknowns, thereby decreasing the compu- 
tational resource demand of the problem. 

The radio frequency (RF) coils within the MRI probe experience an elec- 
tromotive force due to the rotating magnetization. This results in an electric 
current signal which is fed as output to a screen. The signal decays with time 
according to the T2 time constant, hence the name Free Inductance Decay (FID) 
signal [TOT] . 



2.9 Bloch Equations 

The Bloch equations describe the time evolution of the bulk magnetization 
vector, and are therefore also referred to as the equations of motion of magne- 
tization. They are given by, 

^^.(M(.)xB(.)).-^, (23) 



and 



aj -7(M(i)xB(i)),-^^, (24) 



where 7 is the gyromagnetic ratio of the nucleon and B(i) = Bq + Bi(i) 
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2.10 Intravenous MR Contrast Agents 



There are several contrast agents in current clinical use, most of which are 
based either on low molecular weight chelates of gadolinium ion (Gd'^"'") [23lfT6] 
such as Gadolinium-DPTA, or are based on iron oxide (FeO) such as Fe304, 
Fe203, and 7-Fe203 (or gamma phase maghemite) [T71 21 [SHI US]. Based on 
their magnetic properties, they confer contrast of vasculature from background, 
and by proxy of vascularization density, of organs from surround. Specifically, 
the Gd'^"'"-based agents predominantly shorten Ti time, while the FeO-based 
agents predominantly shorten T2 time. Nanoparticle formulations of contrast 
agents are in various stages of development and can be based on gadolinium 
or iron-oxide, or on other agents such as cobalt, nickel, manganese, or cop- 
per ions |109l H] . Optimization of the efficacy of MR contrast agents is an 
active area of research [211 [50l [53l IMl ESI ESI [sl] . The growing array of MR 
contrast applications has included efforts to target tissue or tumor-specific recep- 
tors [MllMlllllinilMllISIinSj- Serious adverse reactions to intravenous MRI 
contrast agents are rare. However in the mid 2000s gadolinium chelates were in- 
creasingly hypothesized as playing a role in nephrogenic systemic fibrosis (NSF) , 
a disease entity unknown prior to 1997. It has since been essentially established 
that Gadodiamide, a gadolinium-containing agent, is associated with an in- 
creased risk of developing NSF in susceptible individuals [TH |S31 jWl 11061 IBB] . 
Superparamagnetic iron oxides have been promoted by some as possible alter- 
natives to gadolinium-based agents in patients at risk for NSF |81) . 



2.11 Gradient-based 3D Spatial Localization 

For 3D localization, a gradient coil generates a magnetic field gradient in the 
bore. Then in accordance with the larmor frequency relation w = jB{r,t), it 
follows that the resonance frequency, w = uj(r,t) acquires space dependence, 
and thereby serves to effectively label the tissue in space. The process of space- 
specific excitation of the sample is called slice selection. The excitation pulse 
frequency is simply chosen as a bandwidth frequency fianking the region-of- 
interest or the slice. The process of scanning is an iteration, however simple 
or complex, through a sequence of slices. And slice thickness is proportional to 
excitation pulse bandwidth. The slice thickness can be decreased by decreasing 
the bandwidth, or alternatively by increasing the magnetic field gradient. For 
instance, given a constant magnetic field gradient, the slice thickness between 
two points ri and r2 is given by, 

-r,-r, (26) 



7[B(r2)-S(ri)] 

Both the bandwidth and magnetic field gradient are Graphical User Interface 
(GUI)- adjustable parameters on current day MRI machines. This gradient- 
based localization technique was Paul Lauterbur's contribution to the develop- 
ment of MRI EE- 



2.12 MRI Machinery 

The MRI is a simple machine in principle. It was was invented by Raymond V. 
Damadian by 1971 [5111 15^ ISTl 1111) . It consists of a magnet, a radio frequency 
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coil, an empty space or bore, and a computer processor. The magnet is for 
applying the static and pulsed magnetic fields. Current day machines typically 
use superconducting electromagnets. These are wire coils cooled to low temper- 
atures using liquid helium and liquid nitrogen |29l |63l [H SH HH [Ml EH ITll] , 
allowing for greatly decreased electrical resistance, and consequently high cur- 
rent loops generating strong magnetic fields. The RF coil typically serves both 
as the generator of the pulsed magnetic field, and as the detector of the FID 
signal. Such duality of function is possible simply because the Maxwell-Faraday 
equation is true in both directions. Several machines also have a separate gra- 
dient coil. The gradient coil has x, y, and z components indicating the spatial 
direction of the gradient field, and together conferring 3D localization. The com- 
puter processor runs the software instructions for pulse sequences and signals 
processing, including image reconstruction. 

2.12.1 Superconductivity 

MRI requires strong magnetic fields. Though the Earth's magnetic field is 
variable with time and is asymmetric between northern and southern hemi- 
spheres [SS], it is estimated to range between 25,000 and 65,000 nanoteslas (0.25 
to 0.65 Gauss) at the surface, with a commonly cited average of 46,000 nanotesla 
(0.46 Gauss) |113j . Current clinical use MRIs commonly range from 1.5 Tesla 
to 7 Teslas. Therefore MRIs may require magnetic fields up to 280,000 times 
the earth's magnetic field. This has been achieved by a shift from permanent 
magnets to superconducting electromagnets, and presents another interdisci- 
plinary link to both theory and practice. The dominant physical theory of su- 
perconductivity is called BCS theory [71 [25l [3 [9] after its formulators, Bardeen, 
Cooper, and Schrieffer, and is based on electron pairing via phonon-exchange 
to form so called Cooper pairs. Though electrons are fermions and subject 
to the Fermi exclusion principle, precluding any two electrons occupying the 
same state. Cooper pairs act like bosons because their pairing and spin summa- 
tion results in integer-spin character, a type of boson-like state. Upon cooling, 
these boson-like pairs can condense, analogous to Bose-Einstein condensation 
and representing a quantum phase transition. Quantum phase transitions have 
historically been challenging to study in a controlled way in the laboratory, but 
have recently been demonstrated by H. T. Mebrahtu et. al. in tunable quantum 
tunnelling studies of Luttinger liquids [701 . In superconductivity, upon transi- 
tion into the condensed state, the coherence of the system imposes a high energy 
penalty on resistance, effectively resulting in a zero resistance state. This is for- 
mally encoded in the theory by an energy gap at low (subcritical) temperatures. 
In superconducting electromagnets, the coil cooling is done through liquid he- 
lium and liquid nitrogen and presents a connection to the chemistry of phase 
transitions, condensed matter and solid state physics, as well as the engineering 
aspects of attaining and maintaining such low temperatures. 

3 MRI Assessment of Acute Ischemic Stroke 

In this section, we review the MRI modalities used in the assessment of acute 
ischemic stroke. 
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Figure 4: ADC map showing hypointensity in the left periventricular MCA 
distributions, consistent with acute ischemic stroke. Note the marked hyperin- 
tensity of the vitreous and cerebrospinal fluid. Both are fluid-filled cavities and 
have a much higher diffusion coefficient than tissue. 

3.1 Diffusion- Weighted Imaging 

The diffusion-weighted MRI (DWl) modality exploits the Brownian motion of 
water molecules within a sample, and the relative phase shift in moving wa- 
ter versus stationary water. The Apparent Diffusion Coefficient (ADC) is a 
parameter which can be mapped to provide diagnostic information. In acute 
ischemic stroke, cytotoxic cellular injury results in axonal edema and a sub- 
sequent decrease in Brownian motion. This is manifested as a hyperintensity 
on DWI and a hypointensity on the corresponding ADC map. Figure ([6]) is a 
DWI showing hyperintensity in the left MCA and PCA distributions, consistent 
with acute ischemic stroke. Figure Q is an ADC map showing hypointensity 
in the left MCA distribution, consistent with acute ischemic stroke. Figure ([s]) 
shows an ADC map with a left parieto-occipital hypointensity reflecting left 
posterior circulation acute ischemic stroke. In the subacute setting, the ADC 
may normalize and even increase due likely to ischemia-associated remodeling 
and loss of structural integrity. The differential for hyperintensity on DWI in- 
cludes hemorrhagic stroke, traumatic brain injury, multiple sclerosis, and brain 
abscesses [HH [TTl [^3 DWI has long been shown in animal models to 

be more efficacious than T2-weighted imaging for the early detection of tran- 
sient cerebral ischemia [73j [77l EH] • Additionally, DWI has been shown to be 
highly efficacious in the early detection of acute subcortical infarctions [10^. 
DWI in conjunction with echo-planar imaging has been shown to effectively dis- 
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Figure 5: ADC map showing hypointcnsity in the PCA distribution, eonsistent 
with a parieto-occipital acute ischemic stroke 

criminate between high grade (high cellularity) and low grade (low cellularity) 
gliomas |104j . Of note, in spite of the relatively high specificity and sensitivity 
of DWI in the early detection of acute ischemic stroke, there is a small subset of 
stroke patients who evade DWI detection in spite of clinically evident stroke-like 
neurological deficits [5] . DWI tractography or diffusion tensor imaging is a form 
of principal component analysis in which the dominant eigendirection is used to 
determine the path of an axonal tract in a given voxel. This method has shown 
potential for further elucidation of neuronal pathways in the brain. 

3.2 Perfusion- Weighted Imaging 

Perfusion-weighted imaging (PWI) is an ordinary differential equations com- 
partment model of blood perfusion through organs. An arterial input function 
(AIF) is determined as input into the particular chosen model. The AIF is 
an impulse, and the model is essentially described by an impulse response or 
Green's function. The computed output are perfusion parameters such as mean 
transit time, time to peak, cerebral blood flow rate, and cerebral blood vol- 
ume. There are various protocols for the perfusion conditions. For instance 
the dynamic susceptibility contrast imaging method uses gadolinium contrast 
to gather local changes in T2* signal in surrounding tissue. T2* is a type of T2 in 
which static dephasing effects are not explicitly RF-canceled, therefore dephas- 
ing from magnetic field inhomogeneities and susceptibility effects contribute to 
a tighter FID envelope (more rapid decay; T2* < T2) than is seen in T^- Other 
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Figure 6: DWI showing hyperintensity in the left PCA and MCA distributions, 
consistent with acute ischemic stroke. 

PWI protocols include arterial spin labeling and blood oxygen level dependent 
labeling. Not surprisingly, PWI results vary with the choice of perfusion model 
and computational method with which it is implemented 11051 [5^ [H^ . 

3.3 Combined Diffusion and Perfusion- Weighted Imaging 

Both PWI and DWI are efficacious in the early detection of cerebral ischemia 
and correlate well with various stroke quantitation scales |108l IHl [S3] . Further- 
more, the combination of diffusion and perfusion-weighted imaging allows for 
assessment of the ischemic penumbra, that watershed region at the boundary of 
infarcted and non-infarcted tissue [80 l [99l fTTSl l98] . It represents tissue which has 
suffered some degree of ischemia, but remains viable and possibly salvageable 
by expedient intervention with thrombolytic therapy [691 1881 118j , hypothermia, 
blood pressure elevation, or other experimental methods under study. Outside 
the penumbra, thrombosis or autoregulation of circulation shuts down perfu- 
sion to infarcted tissue, resulting in a drastic decrease in perfusion. Diffusion is 
simultaneously decreased due to tissue ischemia, hence no PWI-DWI mismatch 
occurs. Within the penumbra however, the tissue is hypoperfused, yet since 
still viable, it has yet to sustain sufficient cytotoxic axonal damage to manifest 
a significant decrease in diffusion. The registration of both imaging modali- 
ties therefore demonstrates a diffusion-perfusion mismatch at the penumbra. 
In set-theoretic terms, the brain territory A with compromised diffusion is a 
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proper subset of the territory S' with compromised perfusion, and the relative 
complement of Z\ in if" is called the ischemic penumbra A, 

A = A"" n'F = ^\A (27) 

One caveat to PWI-DWI mismatch assessment is that when done subjec- 
tively by the human eye, it can be demonstrably unreliable |26) . Hence devel- 
opment of quantitative metrics and learning algorithms are needed in this area. 
For instance, generalized linear models have been used and shown promise |115j . 




Figure 7: 3D time-of-flight image of the carotid vasculature. Axial section of 
neck shown. 

3.4 Magnetic Resonance Spectroscopy 

Magnetic resonance spectroscopy (MRS) is being used in attempts to reliably 
measure the magnetic resonance of metabolites whose concentrations change 
in the setting of acute ischemic stroke. Lactate (LAC) levels are well known 
to increase in ischemia, and is therefore a candidate in development, while N- 
acetyl aspartate (NAA) levels have been shown to decrease in acute stroke, and 
are used in proton spectroscopy studies. Most studies confirming the ischemia- 
associated elevation in LAC and decrease in NAA have been done in the setting 
of hypoxic-ischemic encephalopathy in newborns [2Ql [TOl [67] . In addition to 
lactate, a-Glx and glycine have also been shown to be increased in the as- 
phyxiated neonate [53]. Phosphorus magnetic resonance spectroscopy detects 
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Figure 8: 2D time-of-flight fast spoiled gradient echo sequence (FSPGR) image 
of the carotids. 

changes in the resonance spectrum of energy metaboUtes, and has also shown 
prognostic significance in hypoxic-ischemic brain injury [3ll90). In their current 
form, MRS methods are less sensitive and practical for the assessment of acute 
ischemic stroke than the other magnetic resonance modalities discussed here. 
However, MRS is rapidly finding a place in the routine clinical assessment of 
hypoxic-ischemic encephalopathy in the asphyxiated neonate. Some neonatal 
intensive care units, for instance, now conduct MRS studies on all neonates 
below a certain threshold weight. 

3.5 Blood Oxygen Level-Dependent (BOLD) MRI 

The blood oxygen level-dependent (BOLD) MRI is a magnetic resonance imag- 
ing method that derives contrast from the difference in magnetic properties 
of oxygenated versus deoxygenated blood. Hemoglobin is the molecule which 
carries oxygen in the blood, and is located in red blood cells. The oxygen- 
bound form of hemoglobin is called oxyhemoglobin, while the oxygen-free form 
is called deoxy hemoglobin. Deoxyhemoglobin is paramagnetic and in the pres- 
ence of an applied magnetic field assumes a relatively higher magnetic dipole 
moment than oxyhemoglobin which is a diamagnetic molecule. Deoxygenated 
blood has a higher concentration of deoxyhemoglobin than oxyhemoglobin, and 
this difference is reflected in the MRI signal. r2*-based BOLD MRI has shown 
promise in localizing the penumbra in acute ischemic stroke [HI |3S] . It is also 
used in functional MRI (fMRI) , which is based on the principle that active brain 
areas have higher resource (e.g. oxygen, glucose) demands and higher waste out- 
put 03]. In this context, BOLD has been used to study the brain's behavior 
during sensorimotor recovery following acute ischemic stroke. Specifically, the 
coupling between BOLD and electrical neuroactivity has shed some light on 
the still poorly understood process of spontaneous motor recovery following a 
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Figure 9: MRA showing axial view of the internal carotids and Circle-of- Willis 



stroke [ISl [H [SlIMl 123 • BOLD MRI results can be affected by baseline cir- 
culatory status, and therefore in the research setting near-infrared spectroscopy 
(NIRS) can be used as a control, or as an alternative modality in the clinic 
[1021 [78] . BOLD MRI has been used in assessing the brain's response to hyper- 
capnia [1101 [5] . Hypercapnia can cause changes in multiple variables in the brain 
such as cerebral blood flow, oxygen consumption rate, cerebral blood volume, 
arterial oxygen concentration, and red blood cell volume fraction. However, 
hypercapnia-associated cerebrovascular reactivity has been strongly correlated 
with arterial spin labeling and other PWI surrogates, suggesting that the brain's 
reaction to hypercapnia is dominated by changes in cerebral blood flow j66| . 



3.6 Magnetic Resonance Angiography 

MRA is a set of magnetic resonance-based techniques for imaging the circulatory 
system. Figures ^ and (10) show MRA images of the Circle-of- Willis. MRA 



techniques can be broadly categorized into flow-dependent and flow-independent 
groups. The flow-dependent methods derive contrast from the motion of blood 
in vasculature relative to the static state of surrounding tissue 1331 . Two cur- 
rently well-known and used examples of flow dependent methods are: (i) Phase 
contrast MRA, (PC MRA) [Ml [371 [HJ and (ii) Time-of-Flight MRA (TOF 
MRA). The TOF MRA images can be acquired in either two dimensional (2D 
TOF) or three dimensional (3D TOF) formats. Figure ^ shows a 2D TOF fast 
spoiled gradient echo sequence (FSGR) image of the carotids, and Figure ^ 
shows a an axial section of a 3D TOF image of the carotids. PC MRA exploits 
differences in spin phase of moving blood relative to static surrounding tissue, 
while TOF MRA exploits the difference in excitation pulse {Bi) exposure of 
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Figure 10: MRA lateral view showing internal carotids ascend into the Circle-of- 
Willis and give rise to the middle cerebral arteries. Outline of the ear is visible 
posterior to the internal carotids. Branches of the external carotid arteries are 
also seen as the facial artery anteriorly, and the arteries of the scalp posteriorly. 



flowing blood relative to static surrounding tissue. This difference occurs be- 
cause flowing blood spends less time in the field of exposure, and as a result is 
less spin-saturated then the surrounding tissue. This decreased spin-saturation 
translates into higher intensity signals on spin-echo sequences. 

Flow-independent methods exploit inherent differences in magnetic proper- 
ties of blood relative to surrounding tissue. For instance, fresh blood imaging 
is a method that exploits the longer T2 time constant in blood relative to sur- 
round. This method finds specific utility in cardiac and cardio-cerebral imaging 
by using fast spin echo sequences which can exploit the spin saturation differ- 
ences between systole and diastole [H]. Other examples include susceptibility 
weighted imaging (SWI) and four dimensional dynamic MRA (4D MRA). SWl 
derives contrast from magnetic susceptibility differences between blood and sur- 
round [46: , while 4D MRA uses time-dependent bit-mask subtraction after injec- 
tion of Gadolinium-DPTA or some other pharmacological contrast agent |117j . 
Figures ( 12 ) and ( 13 ) show SWl consistent with acute ischemic infarction in the 
left PGA and MCA distributions. Of note, in the sense that 4D MRA exploits 
the time interval between injection and initial image acquisition, it is arguably 
more flow-dependent than the other methods mentioned here in that category. 



4 The Hydrogen Atom 

Hydrogen is the most commonly imaged nucleus in magnetic resonance imag- 
ing, due largely to the great abundance of water in biological tissue, and to 
the high gyromagnetic ratio of hydrogen. Furthermore, hydrogen is the only 
atom whose schrodinger eigenvalue problem has been exactly solved. Other hy- 
drogenic atoms involve screened potentials and require iterative approximation 
eigenvalue problem solvers such as the Hartree-Fock scheme and its variants. In 
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Figure 11: MRA of the neck showing the carotid vasculature. 



what follows we review the electronic orbital configuration of hydrogen. 

4.1 Electronic Orbital Configuration 

In this section we focus on the electron. Specifically that lone electron in the 
orbit of the hydrogen atom. Wc make this choice because its simplicity allows 
us review the attributes of quantum orbital angular momentum and spin, in a 
context which is not only real, but also highly relevant to magnetic resonance 
imaging, which as noted targets the hydrogen atom. The electron in -'^H has both 
position and angular momentum. The Hamiltonian commutes with the squared 
orbital angular momentum operator, therefore angular momentum eigenstates 
are also energy eigenstates. Of note, the position and angular momentum ob- 
servables do not commute. This is equivalent to Heisenberg's uncertainty princi- 
ple. Hence definite energy or orbital angular momentum states, are represented 
in the position basis as probability density functions, which are synonymous 
with electron clouds or atomic orbitals. 

Orbital angular momentum is the momentum possessed by a body by virtue 
of circular motion, as along an orbit. In quantum mechanics, the square of 
the orbital momentum is a quantized quantity, which can take on only certain 
discrete "allowed" values. This reality and discreteness of the orbital momen- 
tum eigenvahics is a manifestation of the spectral theorem of normal bounded 
operators; given the hermiticity and boundedness of the quantum orbital angu- 
lar momentum operator over a specified finite volume such as the radius of an 
atom; which is itself a manifestation of the negative energy of the electron in 
orbit yielding a so called bound state. The orbital momentum states can be indi- 
rectly represented by the spherical harmonic functions to be derived below. The 
spherical harmonics are the cigenfunctions of the squared quantum orbital angu- 
lar momentum operator in spherical coordinates. Their product with the radial 
wave functions yield the wave function, tjj, of the hydrogen's electron. Where 
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Figure 12: SWI showing signal in the left MCA distribution consistent with 
acute ischemic infarct. Note the left periventricular region of hyperintensity. 
Its hypointense center is consistent with the ischemic core. 



in accordance with the principal postulate of quantum mechanics, ip'^{r, 6,(j)) is 
the probability of finding the electron at a given point (r, 9, (j)). 

In the time-independent case, the Hamiltonian defines the Schrodinger equa- 
tion as follows. 



Hip = Ell;. (28) 

The orbital angular momentum states are the angular portion of the solu- 
tion, V', of the schrodinger equation in spherical coordinates. Substituting the 
expression for the Hamiltonian of a single non-relativistic particle into Equa- 



tion (28) yields, 

- I^V'V-W + V{v)i;{v) = EiPir), (29) 
where r e 3?"^. In spherical coordinates the above equation becomes. 



2m sin t 



dr 



dr 



d_ 
86 



sin( 



86 



1 92 



smf 

V{r)ijj{r) 



EiPir), 



(30) 



where m denotes mass. Invoking the method of separation of variables, we 
assume existence of a solution of the form %p{r) = R{r)Y [6 , cf)) . We illustrate 
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Figure 13: SWI showing signal in the left PCA and MCA distributions consis- 
tent with acute ischemic infarct. Note the bull's eye pattern of hyperintensity 
surrounding some central foci of hypointensity in the occipital lesion. 

the method on the Laplace equation which can be interpreted as a homogeneous 
form of the schrodinger equation, i.e. one for which V{r) = E = 0, 

= 0. (31) 

Multiplying through by r'^/RY yields, 

^ 9 f 2d \ ^, , 1 d f . ..,n 1 

Rd-r\:d-r) ^(^^ ^Y^ede ^ 9^ ) + Y^O W 

= 0, (32) 

which can be separated as, 
and 



21 



We can again assume separability in the form Y{9,(p) = Q{9)^{(j)), and 
substitute into Equation (34 1 above to get, 



19/ d\ 1 d^^ 

sine- m + — = -t (35) 



esinede \ de) ^ ' ^ sm^ e dcj)'^ 

From which we extract the following two separated equations 

1 



^ (36) 



and 



, . 2. sin6i d f . BQ\ 2 
Our homogeneous equation has therefore been separated into three ordinary 



differential equations: Equations (33) for the radial part, Equation (36) for the 
azimuthal part, and Equation (37) for the polar part. Solutions to the radial 
equation are of the form, 

R{1) ^ Ar^ + Br-^'+^^ (38) 

where A and B are constant coefficients, I is a non-negative integer such 
that I > |m|, where m is an integer on the right hand side of the azimuthal 
and polar equations. A regularity constraint at the poles of the sphere yield a 
Sturm-Liouville problem which in turn mandates the form ^ = 1(1 + 1). The 
coefficient A is often set to zero to admit only solutions which vanish at infinity. 
However, this choice is application-specific, and for certain applications it may 
be appropriate to instead set B = 0. 

Solutions to the azimuthal equation are of the form, 

Ce-™'^-HL»e+™'^, (39) 

where C and D are constant coefficients and e is the base of the natural 
logarithm. To obtain solutions to the polar equation, we proceed in a number 



of steps. First we substitute cos 6 x, and recast Equation (37) into, 

sin^ 0-^ + sin 6* cos 6*— + ?(/ + 1)6 sin^ 9 - m^Q = 0, (40) 

Next we compute the derivatives of S under the transformation x — cos 9. 
Employing the chain rule yields, 

de dQdx . dQ 

-d9^d^T9=-''''^^ 

and 



d^Q d f . de\ de . ^d dQ 

= j =-''''^d^~''''^d9d^ 



^sm^9-— -cos9—. (42) 
dx^ dx 

A substitution of the derivatives into the polar equation yields. 
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sm 6 sm — cos —— + sm cos a — sm — — 
\ dx-^ dx J \ dx 

+l{l + l)Qs\T?e-m^Q^Q. (43) 

Dividing through by sin^ and changing variables via cos 6 i-^ x and Q i-^ y, 
yields the associated {m ^ 0) Legendre differential equation, 

whose solutions are given by 

p-^x) = ~ x^y^'^^j^ix^ - 1)'. (45) 

The solution F(^, 9) = ^{(f>)Q{d) to angular portion of the Laplace equation 
is therefore of the form, 

Yi,,n{0, (t>) = iVPr(cos0)e™'^, (46) 
where N is a normalization factor given by, 

N=,^^^±M^ (47) 

Y 4^ (/ + 1)!' ^^'> 

and enabling, 

' Y*„,{n)Yi^„,{n) = 1, (48) 

where = (0, 9) is solid angle. Y; ,„ are called the spherical harmonic 



functions and are discussed some more in the following subsection. Figures ( 14 ) 



to (23 1 are plots of a sampling of spherical harmonics with 1 < ^ < 5. 

For the electronic configuration of the hydrogen atom, the electron experi- 
ences a potential V(r) due to the proton. V(r) is the coulomb potential given 

by, 

Vir)^~-^ (49) 
Substituting this into the generic schrodinger equation yields, 

- |^VV(0 - j^i^ir) = Ei;ir), (50) 
which in spherical coordinates is. 



2m sin 9 



&in9— r^— + — sme* 



dr \ dr J 89 \ 89 J sin 9 8(p 

2 



■0(r) 



47reor 
where m is given by, 



^{r)=E^P{r), (51) 
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Figure 14: 1=1, m=l 



mpni^ 



(52) 



and is the two-body reduced mass mp of the proton and of the electron. 



Equation (51 ) above is an eigenvalue problem which after some rearranging, 



we can solve using the same separation of variables method illustrated above. 
We recast as. 



d_ 
dr 



d \ 2m f re^ 



1 



dr 
d 



sm0d0 



sm 



Er'' 



del sin^ e dS^ 



R{r)Y{e,^) 

R{r)Y{e,^) 
= 



Dividing through by R{r)Y{9, (p) and splitting the operator, we get. 



Y{n) 



1 d 



de 



1 



sin^ 9 



(53) 



(54) 



and 



1 d 

R{r) dr 



dr 



Rir) 



2m 

If 



re 



Er' 



= /(/ + l). 



(55) 



Equation (54) is exactly the same as Equation (35) which we solved above. 



and whose solutions are the spherical harmonic functions, Yz^m- Equation (55) 



above is isomorphic to a generalized Laguerre differential equation, whose so- 
lutions are related to the associated Laguerre polynomials. As mentioned in 
the case of Laplace's equation above, regularity conditions at the boundary 
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Figure 15: 1=2, m=l 



prescribe the separation constant 1(1 + 1). To transform Equation (551 into a 
generalized Laguerre equation, we recast into the following form. 



2mre 



-1(1 + 1) 



R(r) = 0. 



then we proceed with the following sequence of substitutions: 



(56) 



R(r) ^ 



y(r) 



2mE /32 

X 



Executing the first substitution above transforms Equation ( 56 1 into, 
(fv(r) \ 2me2 2'mE 1(1 + 1) 



y(r) = 0. 



Next, the second substitution transforms the above into, 



(Py(T_ 



2me^ 



4 



1(1 + 1) 



y(r) = 0. 



And finally, the third substitution transforms the above into. 



^2(Py{x) 
dx^ 



2/3me2 

ATTXenh? 



A 



^""1(1 + 1) 



Dividing through by 0^ yields. 



<Py(x) 

dx-^ 



1 



2me 



4 AirPe^h'^x 



y(x) = 0. 



y(x) = 0. 



(57) 



(58) 



(59) 



(60) 



(61) 
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Figure 16: 1=3, m=l 



Next, we make the substitutions. 



e - 1 



and 



2me 
47r/3eo/i2 



2j + k + l 



(62) 



(63) 



which transform Equation (61 ) into the fohowing generahzed Laguerre equa- 



tion 



1 2j + k 

— h — 

4 2x 



1 _ P-1 



4(3^) =0, 



(64) 



whose solutions are readily verified to be the following associated Laguerre 
functions. 



where L*^ are the associated Laguerre polynomials. 
From Equation (62 1 we see that, 



2/ + 1. 



which we substitute into Equation ( 63 1 to yield, 
2me^ 2j + (2^ + 1) + 1 



j + l + l = n. 



(65) 



(66) 



(67) 



47r;3eoS2 2 

n is the principal (or radial) quantum number. The corresponding energy 
eigenvalues, are obtained by substituting /J^ into and solving for E, we 
get, 
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Figure 17: l=3,m=2 



where ao is the Bohr radius and is given by, 



(68) 



(69) 



The value of £'„ in the ground state (n = 1) is called the Rydberg constant. 
Its value is shown in Table ^ along with that of other physical constants 
pertinent to magnetic resonance imaging. The displayed values are from the 
2010 Committee for Data on Science and Technology (CODATA) recommenda- 
tions [75] , and Ur is the relative standard uncertainty. 



To transform the radial equation solution, Equation (65), into a form in 
terms of radial distance r, and principal and azimuthal quantum numbers n 
and Z, we note the following relations: 



From Equation (661 we get, 



1 



(2/ + 1) 



l + h 



(70) 



from Equation ( 67 1 we get 



3 + l + 



and from Equations (57 1 and (68) we get 



4 

Therefore, 



2mEn 



2m 



2maQn'^ 



X = l3r 



2r 



/3 = 



(71) 



(72) 



(73) 
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Figure 18: l=3,m=3 
Substituting the above derived expressions of x, j, and k into the associated 



Laguerre function, Equation (65), we get. 

Next we make the substitution R ~ y{r)/r, which yields the radial solutions 
of the hydrogen atom, 



\aon / 



\aonJ 



(75) 



where the constant C — Following multiplication by the spherical har- 
monics and normalization, we obtain the exact solution of the time-independent 
Schrodinger equation of the hydrogen atom. 



{r,9,(p\n,l,m) = 'I'(r, 6*, c 



2 ^'i!i^l^e-^/-r^VL-+,\m 



nao / 2rt[(n -I- 1)!J'^ \nao J \nao J 

where the principal, azimuthal, and magnetic quantum numbers (n, I, m) 
take the following values. 
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Figure 19: 1=4, m=l 



n = 1,2,3... 
; = 0,l,2,...n-l 
m = -I, -I + 1, ...0, ...I - 1,; 

4.2 Spherical Harmonics and Radial Wave Functions 

In this section, we review the spherical harmonics and the radial wave functions. 
As noted above, the spherical harmonics, YJ™' , are characterized by an azimuthal 
and a magnetic quantum number, I and m; respectively. Table Q shows the 
spherical harmonics with I = 0,1,2, and 3. Formalism of state representation by 
the spherical harmonics is as follows, 

\l,m)=Yi^m (77) 

{0,^\l,m) =Yi^miO,(b) (78) 

where, < < tt is the polar angle and < </) < 27r is the azimuthal 
angle. The Spherical harmonic representing rotation quantum numbers \l,m) 
at position state 16*, (j)), is given by 

^ ^ ^^4^^^ ll^y Pricose)e-t (79) 
4.2.1 Operator Representation 

The spherical harmonics are a complete set of orthonormal functions over the 
unit sphere. In particular, they are eigenfunctions of the square of the or- 
bital angular momentum operator, L^. And are thereby representations of the 
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Figure 20: l=4,m=2 

allowed discrete states of angular momentum. can be arrived at by rep- 
resenting Laplace's equation in spherical coordinates and considering only its 
angular portion. Classically, angular momentum is: 

L = r X p, (80) 

where r is the position vector and p is the momentum vector. 

By analogy, the quantum orbital angular momentum operator is given by, 

\. = -ih{vxV), (81) 

where V is the gradient operator, and we have used the momentum operator 
of quantum mechanics, 

p = -ihV. (82) 
In spherical coordinates, is then represented by, 

and similarly, Lx, Ly, and Lz are represented as follows, 

f d d\ 
Lx ^ ih [sin cot 9 cos cj)—j , (84) 

f d d\ 

Ly — ih [ — cos ch— + cot 9 sin ch— , (85) 
\ o9 ocp J 

Given the expressions above for Yi^ra{9,4'), L^, and L^, it is readily shown 
that, 
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Figure 21: l=4,m=3 



L^\l,m) =h^l{l + l)\l,m), (87) 

and 

L^|/,m) = nm\l,m). (88) 

The square of the radial wave functions is the probabiUty the electron is 
located a distance r from the nucleus. Table ^ shows the radial wave functions 
for n — 1,2, and 3. 

4.2.2 Commutation Relations and Ladder Operators 

[Li,Lj] = iheijkLk, (89) 
where Cijk is the Levi-Civita symbol given by, 

{+1 if (i,j,k) is an even permutation of (1,2,3), 
— 1 if (i,j,k) is an odd permutation of (1,2,3), (90) 
if any index is repeated, 

and i,j, k can have values of x, y, or z. 

[L\Lj]=Q fori = a;,y,z (91) 

The ladder operators act to increase or decrease the m quantum number of 
a state and are given by, 
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Figure 22: 1=4, m=4 



'V 



(92) 



It is readily verified that 



L±\l,m) = h^/{l ^m){l + l±m)\l,m± 1), 



(93) 



[L,,L±]^±hL±, 



(94) 



and 



(95) 



5 Intrinsic Spin 

Spin is an intrinsic quantum property of elementary particles such as the elec- 
tron and the quark. Composite particles such as the neutron, proton, and even 
atoms and molecules also possess spin by virtue of their composition from their 
elementary particle constituents. The term spin is itself a misnomer, as a phys- 
ically spinning object about an axis is not sufficient to account for the observed 
magnetic moments. The mathematics of spin was worked out by Wolfgang Pauli, 
who either astutely or serendipitously neglected to name it. This was wise, as 
later insight elucidated spin as an intrinsic quantum mechanical property with 
no classical correlate. 
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Figure 23: l=5,m=2 

Spin follows essentially the same mathematics as outlined above for orbital 
angular momentum. This is by virtue of the isomorphism of the respective Lie 
Algebras of S0(3) and SU(2) groups. And is elaborated further in section ([t]) 
below. One notable difference is that half-integer eigenvalues are allowed for 
spin, while orbital angular momentum admits only integer values. The spin 
algebra is encapsulated in the commutation relations as follows, 

[Si, S-j] = iheijkSk, (96) 

[S^,S,]^0 for J = a:, 2/, z, (97) 

S±=S^±iSy, (98) 

S±\s,m) = h^/{sTm){s + l±m)\s,m±l), (99) 

[S,,S±]=±hS±, (100) 

and 

[S+,S-]=2hS^. (101) 

6 Addition and the Clebsch-Gordan Coefficients 

The addition of spin and or of orbital angular momentum refers to the process 
of determining the net spin and or orbital angular momentum of a system of 
particles. For example, the spin of a hydrogen atom in its ground state, i.e. the 
I = state with zero orbital angular momentum, constitutes the composite spin 
of the electron and the proton. Similarly, the nuclear spin, /, is the composition 
of the spins of protons and neutrons in the nucleus of an atom. 
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Table 2: Physical constants of magnetic resonance imaging 



Quantity 


Symbol 


Value 


Units 




vacuum light speed 


C, Co 


299792458 


m s~^ 


exact 






Planctc constant 


h 


6.62606957(29) x 10"^* 


J s 


4.4 X 


10" 


8 


Reduced Planck 


h/2TT, h 


1.054571726(47) x 10"'*" 


J s 


4.4 X 


10" 


8 


Rydberg constant 




10973731.568539(55) 


m-i 


5.0 X 


10" 


12 


Rydberg energy 


Ry, hcR^ 


13.60569253(30) 


eV 


2.0 X 


10" 


8 


Bohr radius 


ao 


0.52917721092(17) x 10""-" 


m 


3.2 X 


10" 


10 


Bohr magneton 


fJ-B 


927.400968(20) x 10"^'' 


J T-^ 


2.2 X 


10" 


8 


Nuclear magneton 




5.05078353(11) x 10"^^ 


J T-^ 


2.2 X 


10" 


8 


Electron g-factor 


9e 


-2.00231930436153(53) 




2.6 X 


10" 


13 


Proton g-factor 


Qp 


5.585694713(46) 




8.2 X 


10" 


9 


e gyromagnetic ratio 


Je 


1.760859708(39) X lO" 


s-^T-i 


2.2 X 


10" 


8 






28024.95266(62) 


MHz T-^ 


2.2 X 


10" 


8 


p gyromagnetic ratio 


7p 


2.675222005(63) X 10* 


s-^T-i 


2.4 X 


10" 


8 




7p/27r 


42.5774806(10) 


MHz T-^ 


2.4 X 


10" 


8 


e magnetic moment 


Me 


-928.476430(21) x lO"^" 


J 


2.2 X 


10" 


8 


p magnetic moment 


Mp 


1.410606743(33) x lO"^" 


J T-^ 


2.4 X 


10" 


8 


Elementary charge 


e 


1.602176565(35) X 10"^^ 


C 


2.2 X 


10" 


8 


Electric constant 


«o 


8.854187817... X IQ-^^ 


F m-i 


exact 






Electron mass 


me 


9.10938291(40) X 10"^^ 


Kg 


4.4 X 


10" 


8 


Proton mass 


mp 


1.672621777(74) x lO"^'' 


Kg 


4.4 X 


10" 


8 


Neutron mass 


rrin 


1.674927351(74) x lO"^'' 


Kg 


4.4 X 


10" 


8 


mp to me ratio 


mp/me 


1836.15267245(75) 


Kg 


4.1 X 


10" 


10 


Avogadro constant 


Na 


6.02214129(27) X 10"^^ 


mor^ 


4.4 X 


10" 


8 


Molar gas const 


R. 


8.3144621(75) 


J mol-^ 


4.4 X 


10" 


8 


Boltzmann constant 


R/Na, k 


1.3806488(13) X lO"^'* 


J K-i 


9.1 X 


10" 


7 


Electron volt 


(J/C)e, eV 


1.602176565(35) X 10"^^ 


J 


2.2 X 


10" 


8 



Mathematically, the process is described by a change of basis. A change 
from the product space basis ( ji , mi | {j2 , m2 1 to the "net sum" space basis {j, m\ , 
where m = mi +m2 and \ji — J2I < j ji + J2- The representation of the sum 
space eigenstates in terms of the product space eigenstates are referred to as 
the Clebsch-Gordan coefficients. 

31 32 

\j,m)= ^ ^ C{mi,m2,m)\ji,mi)\j2,m2), (102) 

mi=-ii m2=-j2 

where C{mi,m2,m) are the Clebsch-Gordan coefficients and are given by, 

C(mi,m2,m) = (j, m|jij2, mim2), (103) 
and we have notationally represented the uncoupled product space basis by, 

\ji,mi)\j2,m2) ■■= \jij2,mim2). (104) 
6.1 Addition Algorithm 

Here we review a procedural description of the addition of angular momenta. 
Consider a hydrogen atom in the ground state. The spin contributors are 
the lone electron in the 1=0 shell, and a proton which is the sole constituent of 
the nucleus. The possible configurations of spin are: 

in), I U), I ;t), and I tt) (105) 

where f indicates spin +1/2 (spin up) and I indicates spin —1/2 (spin down). 
And where each of the above states are prepared by measuring Sf, S2, Sz,i 
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Table 3: Spherical Harmonics for I = 0, 1, 2, and 3. 
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and S'2_2 for each state. We are interested in the combined state of the two 
particle system, and this requires the quantum mechanical addition of angular 
momenta, or in this example, spin. The net result is a state \s,m) whose spin 
can be determined by measuring 5*^ and whose z-component of spin can be 
determined by measuring 5*^, where S = Si + S2 and 5^ — Sz^i + 8^,2- It 
follows that m = mi +m2, and therefore the possible values of m are 1, 0, 0, -1. 
This suggests a triplet state (s=l) and a singlet state (s=0). The triplet state 
corresponds to m=l,0,-l, and is given by, 

|i,i> ^ in> 

|i,o) ^ ;^(in) + ut)) (106) 

and the singlet state corresponds to, 

|0,0) .^-^(|n>-Ut», (107) 
where upon designation of the top assignment |1, 1), each of the other states 



are obtained by application of the ladder operators defined in Equation ( 99 1 . 
For example, 

S'_|l,l) =fi\/2|l,0) (108) 
To change to a representation in the uncoupled basis, we write. 
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Table 4: Radial Wave Functions for n = 0, 1, 2, and 3. 
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and therefore, 



which yields, 



ii,o) = ^(5i,_in> + 52,-in», 



|i,o> = ^(Ut) + in», 



(109) 



(110) 



(111) 



Finally to obtain the singlet state, we simply orthogonalize the above, yield- 



ing, 



|0,0) 



hV2 



(I ;t) - 1 n)) 



(112) 



The above procedure is easily carried out by a computer for any combination 
of particle spins. For many particle systems, for instance in the computation of 
nuclear spins or non-ground state configurations, the associativity property 
is used and the above description applies. In the above example, the factors 
are the Clebsch-Gordan coefficients. There are several openly available imple- 
mentations of Clebsch-Gordan coefficient calculators, in addition to tabulations 
in handbooks of physics formulae [1141 [T] . 



7 Group Theory: SO(3), SU(2), and SU(3) 

Orbital angular momentum and spin are the source of magnetic resonance. Their 
abstract mathematical description is group theoretic, and extends naturally into 
other fundamental physics such as the strong interaction of quarks in nucle- 
ons and nucleons in nuclei. The Frenchman Henri Poincare commented that 
"mathematicians do not study objects, but the relationships between objects" . 
Relationships between symmetry groups have indeed been the essential device 
for probing the unseen in the realm of high energy particle physics. 

For spin 1 /2 particles such as the electron and the quark, the spin observable 
can be represented by the SU(2) group. SU(N) is the group of N x N unitary 
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matrices with unit determinant, in which the group operation is matrix- matrix 
multiphcation. Similarly SO(N) is the group of N x N orthogonal matrices of 
unit determinant, with matrix multiplication as group operator. Analogous to 
the SU(2) representation of spin, quantum orbital angular momentum can be 
represented by the S0(3) group. Both groups share a Lie Algebra, encoded in 
the commutation relations shown above, because they are isomorphic to each 
other. In particular, SU(2) is a double cover of S0(3). The SU(3) group models 
the interactors of the strong force in Quantum Chromodynamics (QCD). 

7.1 SO(3) 

A counterclockwise rotation by an angle (f) about the x, y, or z-axes, can respec- 
tively be represented by 



cos sin ^ 

RyW =1 1 I , (113) 

— sin cos (j) 

cos (t> — sin 
Rz{4') ^ I sin(^ cos(/) 
1 

The corresponding infinitesimal generators, G^, Gy, G^ respectively are: 



(114) 

The corresponding Lie Algebra can then be readily shown to be: 

[Gi, Gj] = ihcijkGk (115) 

7.2 SU(2) 

SU{2) = \( I ~f ) : a,6e C,a2 + 62 = l| (116) 






b a 

It follows that the generators of SU(2) are the Pauli Matrices given by, 



^1 = ( J J ) , ^2 = ( ° ) , ^3 = ( J j , (117) 

The SU(2) Lie Algebra is then prescribed by the following commutation and 
anti-commutation relations, 

[o-a,C^6] = 2iea6cO-c, (118) 
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and 

{aa,ab} = 2Sab- 1 + ieabcO-c, (119) 
which combine to give: 

O-aO-b = Sab ■ I + ieabcCTc, (120) 

where 6ab is the Kronecker delta and {x, y} := xy+yx is the anti-commutator. 

7.2.1 SU(2) is Isomorphic to the 3-Sphere 

Given U e SU{2), such that, 

\ b a* 

We define an assignment, (, such that 

^ ja^xo + ix, ^ ^^21) 
I 6 —X2 + ixi 

where a;o,a;i,a;2,a;3 e 9?. It follows that, 

U = xqI + ixpGp, (122) 

where I is the 2x2 identity matrix, Gp are the Pauli matrices, and we have 
used the Einstein summation notation over p. 

det{U) = 1^ + b'^ = 1^ xl + XpXp = 1 ^ x e 5^ (123) 

C is a bijective homomorphism between SU{2) and S^. And this shows 
SU{2) is isomorphic to S^. 

7.3 SU(3) 

The generators of the SU(3) group are given by the Gcll-Mann matrices: 
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(124) 



38 



Then defining Tj = the commutation relations are, 

8 

[T,,n]=iY,fabcTc, (125) 

c=l 

and 

4 ^ 

{Ta, n} = -6ab + 2 ^ dabcTc, (126) 
c=l 

where /obc and date are the antisymmetric and symmetric structure constants 
of SU(3) respectively. 



8 Summary 

Magnetic resonance imaging is an application of quantum mechanics which has 
revolutionized the practice of medicine. Electrons, protons, and neutrons are 
magnets by virtue of their spin and orbital angular momentum. Tissue is made 
of ensembles of such subatomic particles, and can be imaged by sensing their 
responses to applied magnetic fields. As discussed in this paper, a wide array of 
MRI modalities exist, each of which derives contrast by specific perturbations 
of the magnetization vector. The potential of MRI in the diagnosis and civen 
the treatment of various diseases is far from fully realized. A quantum level un- 
derstanding of magnetic resonance technology is essential for further innovation 
in this field. 



39 



Acknowledgement 



The author thanks his friends Henok T. Mebrahtu and Peter Q. Blair for read- 
ing this manuscript and providing helpful feedback. And he thanks them for 
many enjoyable physics conversations over coffee. He thanks Daniel Ennis for 
making available his MATLAB code for plotting spherical harmonics. He thanks 
his lovely wife Lisa M. Odaibo for a helpful conversation on the clinical use of 
magnetic resonance spectroscopy in hypoxic- ischemic injury assessment in the 
neonate. And he especially thanks her for her love and support. The author 
thanks his entire family, and celebrates his father, Dr. Stephen K. Odaibo, on 
his birthday today. He thanks him for being a kind wonderful father and ex- 
ample. He thanks those who contributed and are contributing to his education, 
including: Dr. Robert A. Copeland Jr., Dr. Leslie Jones, Dr. Janine Smith- 
Marshall, Dr. Bilal Khan, Dr. David Katz, Dr. Earl Kidwell, Dr. William 
Deegan, Dr. Brian Brooks, Dr. Natalie Afshari, Dr. Isaac O. Karikari, Dr. 
Xiaobai Sun, Dr. Mark Dewhirst, Dr. Carlo Tomasi, Dr. John Kirkpatrick, Dr. 
Nicole Larrier, Dr. Brenda Armstrong, Dr. Phil Goodman, Dr. Joanne Wilson 
and Dr. Ken Wilson, Dr. Srinivasan Mukunduan Jr., Dr. David Simel, Dr. 
Danny O. Jacobs, Dr. Timothy Elston, Drs. Akwari, Drs. Haynes, Dr. John 
Mayer, Dr. Lex Oversteegen, Dr. Peter V. O'Neil, Dr. Anne Cusic, Dr. Henry 
van den Bedem, Dr. James R. Ward, Dr. Marius Nkashama, Dr. Michael W. 
Quick, Dr. Robert J. Lefkowitz, Dr. Arlie Fetters, and several other excellent 
educators and role models not mentioned here. 



40 



Appendix A 



Signals processing plays an important role in magnetic resonance imaging, and 
is briefly reviewed in this appendix. The free inductance decay signal is often 
collected in the frequency domain and Fourier transformed into the time do- 
main [35] . Also with a gradual shift towards MRI scanners of higher magnetic 
field, the signal-to-noise ratio requires increasingly more effective noise filter- 
ing algorithms. As signal generation and collection methods become increas- 
ingly more sophisticated e.g via spirals and other geometric sequences, non- 
uniform sampling, and exotic Bi pulse sequence configurations, signals process- 
ing methods need to advance in tandem to address problems arising. Computing 
speed also takes on new significance with such developments. Indeed imaging 
and graphics applications have significantly driven the interest and funding for 
fast algorithms and computing. Notable advancements include the fast Fourier 
transform, the fast-multipole method of Greengard and Rokhlin [101 ST] , and 
parallel computing hardware such as the Compute Unified Device Architecture 
(CUDA) multi-core processors by Nvidia for graphics processing j86j. Further- 
more, a new field is taking form at the interface of high performance software 
algorithms and multi-core hardware processors [HSl [Ml [Ml HSl [HZl [ST] . 
We take a brief look at some of the basics of signals processing below, 

A.l Continuous and Discrete Fourier Transforms 

The Fourier transform converts a time domain signal to its corresponding fre- 
quency domain signal, and is given by, 

/oo 
f{t)e'^-'dt, (127) 
-OO 

while the inverse Fourier transform converts a frequency domain signal to 
its corresponding time domain signal, and is given by, 

/oo 
/(c.)e-*d^. (128) 
-oo 

In practice, the discrete versions of the above equations are used instead, 
and give good results provided that signal data is sampled at an appropriate 
frequency, the Nyquist frequency. Given an MRI signal of N sample points, 
/[n], where n = 0, .., iV — 1, the frequency domain representation is given by the 
discrete Fourier transform, 

N-l 

and the discrete inverse Fourier transform is given by, 

1 

fc=0 

The above Fourier formulas readily allow extension to two and three dimen- 
sions as well as various other customizations, tailored to suite the particular 
MRI application and data format. 
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A. 2 MRI Sampling and the Aliasing Problem 




10 



-10 



-10 



Figure 24: 3D sine function 



The Whittaker- Kotelnikov- Shannon- Raabe- Someya- Nyquist theorem pre- 
scribes a lower bound on the sampUng frequency necessary for perfect recon- 
struction of a band-limited signal. It states that given a band limited signal 
g{t) whose frequency domain representation is such that |5(/)| < B for all / 
and some B, then by sampling at a rate fs = 2B, the image can be perfectly 
reconstructed. The sampling interval is T = l//s, and the discrete time signal 
is expressed as (7[?t.T'] where n is an integer. = ujigjl-K and is in units of hertz, 
hence T is in seconds. The perfect reconstruction is given by the Whittaker- 
Shannon interpolation formula, 



In practice, low-pass filtering methods can be used to decrease the amplitude 
of frequency components which exceed the effective bandlimit. 

The aliasing problem arises when the sampling rate /s is less than twice the 
bandlimit B. The effect of this is that the reconstructed image is an alias of 
the actual image. The Discrete Time Fourier Transform (DTFT) of a signal is 
given by. 




n— — OO 



(131) 



where sine is the normalized sampling function given by. 




(132) 



42 



X{w)= x[n]e-'"", (133) 

n— — oo 

and in the case of an under-sampled MRI signal, the DTFT matches that 
of the alias. Therefore the reconstructed image is merely an alias. This phe- 
nomenon can be described more formally by invoking the poisson summation 
formula, 

oo oo 

X,{f) ■.= Y,X{f-kh)=Y,T-x{nT)e-'-'^f^'^, (134) 

— oo — oo 

where the middle term above is the periodic summation. For a bandlim- 
ited under-sampled signal, increasing the sampling rate to satisfy the Nyquist 
criterion will restore perfect reconstruction. 

A.3 The Convolution Theorem 

The convolution theorem facilitates application-specific image processing before 
or after reconstruction in either k-space or time domain. The theorem is given 

by, 

n{f*9)m = f{^)-9{^), (135) 

and states that the Fourier transform of the convolution of two functions 
equals the product of the Fourier transforms of both functions. Where the 
convolution is defined as, 

/oo 
f{T)g{t - T)dT. (136) 
-CSO 
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